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ABSTRACT 



', Many models of Gamma Ray Bursts (GRBs) invoke a central engine consisting of 

I a black hole of a few solar masses accreting matter from a disk at a rate of a fraction 

I to a few solar masses per second. Popham et al. and Narayan et al. have shown 

that, for M > O.1M0S~^, accretion proceeds via neutrino cooling and neutrinos can 
carry away a significant amount of energy from the inner regions of the disks. We 
' improve on these calculations by including a simple prescription for neutrino transfer 

and neutrino opacities in such regions. We find that the flows become optically thick 
. to neutrinos inside a radius i? ~ 6 — 40i?s for M in the range of 0.1 — IOMqs"^, 

^ \ where FL, is the black hole Schwarzchild radius. Most of the neutrino emission comes 

from outside this region and, the neutrino luminosity stays roughly constant at a 



I value Ly ~ lO^^ergs"^ . We show that, for M > IMqs"^, neutrinos are sufficiently 

' trapped that energy advection becomes the dominant cooling mechanism in the flow. 

■ These results imply that vv annihilation in hyperaccreting black holes is an inefficient 

■ mechanism for liberating large amounts of energy. Extraction of rotational energy by 
Q^i magnetic processes remains the most viable mechanism. 

o 

^ ■ Subject headings: accretion, accretion disks — black hole physics — gamma rays: 

bursts — radiation mechanisms: thermal 

^ ■ 

•i-H - 

%■ 

. rz , 1. Introduction 

Most popular models of Gamma Ray Bursts (GRBs) invoke a binary merger or a collapse 
involving compact objects. In particular, these include mergers of double neutron star binaries 
(Eichler et al. 1989; Narayan, Paczynski & Piran 1992; Ruffert &: Janka 1999), mergers of a 
neutron star with a black hole (Paczynski 1991; Narayan et al. 1992; Ruffert & Janka 2001 and 
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references therein), of helium star with a black hole (Fryer &; Woosley 1998), "collapsars" or "failed 
supernovae" (Woosley 1993; Paczynski, 1998; MacFadyen & Woosley 1999; MacFadycn, Woosley 
& Heger 2001) and "supranovae" (Vietri & Stella 1998). All of the above scenarios lead to the 
formation of a black hole with a debris torus or disk around it. (The only exceptions are models 
in which the GRB energy is provided by the magnetic and rotational energy of the newly formed 
neutron star; e.g.; Usov 1992). In order to understand how the extraordinary amount of energy 
characteristic of GRBs can be extracted, we are motivated to further examine the properties of 
such compact and massive disks around black holes. 

Accretion models in the context of GRBs have been recently discussed by Popham, Woosley 
& Fryer (1999; hereafter PWF), Narayan, Piran & Kumar (2001; hereafter NPK) and Kohri & 
Mineshige (2002). The typical mass accretion rates in GRB models are extremely high, of the 
order of a fraction of solar mass up to a few solar masses per second. Under such conditions, 
the gas photon opacities are also very high and radiation becomes trapped (see e.g.; Katz 1977; 
Begelman 1978; Abramowicz et al. 1988). However, at sufficiently high mass accretion rates, 
although energy advection remains important in the outer parts, the disk becomes dense and hot 
enough in the inner regions to cool via neutrino emission. For this reason PWF named these 
disks neutrino-dominated accretion flows (NDAFs). This regime is of particular interest, because 
neutrinos can, in principle (see e.g. NPK; Ruffert &: Janka 1999), tap the thermal energy of the 
disk produced by viscous dissipation and liberate large amounts of its binding energy (via the 
vv e^e" process in regions of low baryon density). For this mechanism to be efficient, though, 
the neutrinos must escape before being advected into the black hole. 

In this paper we investigate the effects of neutrino transport within the context of NDAFs. 
By using a simple prescription to account for neutrino scattering and absorptive opacities, we find 
that, for accretion rates M > 1 Mq s~^ , the gas becomes increasingly more opaque and neutrinos 
become trapped. As a result, energy advection becomes the dominant cooling mechanism in the 
inner regions of the flows. We show that, as M increases, the disk emitting surface moves further 
out in radius and the neutrino luminosity of the flow remains nearly constant. We show that the 
accretion flow luminosity plateaus as it approaches the neutrino Eddington limit of the inner disk 
and as the neutrino cooling efficiency decreases. 

NPK and PWF also noted the importance of neutrino opacity at M > 1 M0S~^ but did 
not allow for neutrino transport effects in their models. Our work therefore complements the 
earlier studies of NDAFs. A detailed treatment of neutrino transfer was included in the numerical 
simulations carried out by Ruffert & Janka (1999), where accretion tori are formed as a result of 
neutron star merging. Consistent with our findings, these authors also show that opacities can be 
high in such tori. However, their results are specific to the parameter space covered by the merger 
model. 

In §2 we identify the dominant neutrino opacity sources and outUne the basic equations we use 
and the approximations we make to the neutrino transfer problem. In §3 we present our numerical 
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results. We delineate the regions of parameter space in accretion rate and radius where the flow 
becomes optically thick to neutrino emission. In §4 we discuss the stability of the accretion flow. 
Finally, in §5, we compute the flow luminosity of our model, compare it with the derived neutrino 
Eddington limit and discuss various implications of our results for GRB energetics. 



2. The physical model 

2.1. Neutrino emission and photodisintegration 

PWF and NPK have calculated models corresponding to steady-state accretion around 
a black hole with a high mass accretion rate. Their work showed that for the full range of 
accretion rates of interest (0.01 — 10Mqs~^ ), and outside a certain radius R ~ lOO-Rg (where 
i?s = 2GM/c^ = 8.85 X 10^ m3 cm is the Schwarzchild radius of a relativistic compact object of 
mass M = 2>m^ Mq), the disk is advection dominated. Inside this radius, and for M > 0.1 Mqs~^, 
the temperature and density of the gas become high enough that neutrino cooling takes over fairly 
abruptly. 

In such accretion flows neutrinos are generated both by neutronization and by thermal 
emission. The most significant thermal processes are: 

(1) Electron-positron pair annihilation (e~ -|- e"*" — > Vi + i>i where i represents both 
electron-type neutrinos, Ve-i^e-, and heavy-lepton neutrinos f^, z^^, i/,-, P,-). The neutrino cooling 
rate per unit volume due to electron-positron pair annihilation is roughly given by (we take the 
PWF approximation to the Itoh et al. 1989; 1990 results): 

where Tn = T/IO^^ K. We assume an equilibrium mixture of e, and r neutrinos and 
antineutrinos^ so that the total cooling rate, when all species are included, is simply obtained by 
multiplying equation (1) by Ni, = 3, the number of neutrino flavors. 

(2) Nucleon-nucleon bremsstrahlung (n + n — > n + n + Vi + Vi). Bremsstrahlung cooling is 
represented by (Hannestad & Raffelt 1998; Kohri & Mineshige 2002): 

Com ^ lO^VfoTfi.S erg cm-3 s"! , (2) 

where piQ is the scaled density, pio = /o/10^° gcm~'^ . 

(3) Plasmon decay. This is the decay rate of the transverse plasmons, which are normal 
photons interacting with the electron gas through 7 — > z/g + i^e, and is estimated by Ruffert, 



^Thermal equilibrium is quickly established (in < 1 ms) between the neutrinos and the matter via weak, neutral 
current annihilation of e^e'^ pairs (e.g.; Bethe, Applegate and Brown 1980; Salpeter & Shapiro 1981). 
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Janka & Schafer (1996) to be 

^^la^mon ^ >< IO^'TiV^ exp-> (1 + 7^) (^2 + erg cm-3 s'^ , (3) 

where jp = 5.565 x 10~^-\/(7r^ + ?7e)/3 and = fJ^e/kT, and we solve for /Xg, the electron chemical 
potential, using equation (24) in Kohri & Mineshige (2002). 

The second type of neutrino cooling is due to neutronization reactions. The most significant of 
these in this context is electron-positron pair capture on nuclei (p+e~ — > n+u^.; n+e+ — > v'c+p), 
also know as the URCA process. The cooling rate per unit volume is given by (see also PWF): 

%N = %-p + %+n - 9-0 X lO^^^pioTfiXnuc erg cm-3 g-i ^ (4) 

where X^uc is the mass fraction of free nucleons approximately given by (e.g.; PWF; Qian Sz 
Woosley 1996): 

Xnuc ~ 34.8p^o'^^^^n ^ exp(-0.61/rn ) , (5) 

with an upper bound of unity. This prescription takes into account the transition from nucleon 

to a-particlcs which effectively shuts off the URCA process. Xnuc is typically very close to zero 
in the outer disk but photodisintegration breaks down a particles in neutrons and protons once 
T reaches about lO^'^ K. The photodisintegration process is also responsible for cooling the gas 
according to (PWF): 

— 2Q _o -1 , , 

9photo-10 P^^""^ erg cm •^s , (6) 
where v is the disk radial velocity defined in §2.3. 

At low optical depths the emission of neutrinos can be computed directly from the rates given 
in equations (1), (2), (3) and (4). NPK and PWF take into account neutrino cooling as described 
in equations (1) and (4) and follow the optically-thin approach in their computations. 



2.2. Neutrino opacities 

Each neutrino emission process, equations (l)-(4), has an inverse process corresponding to 
absorption. In addition, scattering impedes the free escape of neutrinos from the disk. The 
optical depths for neutrinos are therefore given by the inverse of the processes listed above plus a 
contribution from neutral-current scatterings off nucleons. 

The inverse process to equation (1), i.e. the interaction of neutrinos with one another, gives 
rise to an optical depth: 



4(7/8)(7r4 



2.5 X IQ-^T^^H, (7) 
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where the scale height H is defined in §2.3. The term {7/8)aT'^ is the Fermi-Dirac blackbody 
luminosity and is the same for each neutrino flavor, t'l^; a is the Stefan-Boltzmann constant. 

The inverse process to equation (4), absorption onto protons or onto neutrons, leads to a 
neutrino absorptive optical depth in the flow given by (see also NPK), 

ra,eN « 4(^/8^fr4 ^ IQ-^TfiX^.e^loi^ • (8) 

Similarly, the inverse of bremsstrahlung and plasmon processes (given by Eqs. (2) and (3)) 
lead to two additional sources of opacity, Ta,brcmi ''"a.piasmoni which can be worked out in the same 
way as for the previous two. In accordance with NPK and PWF, we find that electron-positron 
pair annihilation and electron pair capture on nuclei are always the dominant emission mechanism. 
Absorption via the inverse of equations (2) and (3) gives rise to much smaller absorptive optical 
depths than from equations (7) and (8). Although all absorption opacity sources are included in 
our code, we will mainly discuss those derived above. 

A far more important opacity source (particularly for heavy-lepton neutrinos) comes from 
neutral-scattering off nucleons {i^i + {n,p} — > + {n,p}). The cross section for momentum 
transfer for neutrino-nucleon scattering is given by 

/ E \ 

CTs.AT = C's,AfO"o ^ , (9) 

where E is the center of mass neutrino energy, trie is the electron rest mass, c the speed of 
light, (To = 1.76 X 10~^^ cm^ , and Cs^n = (1 + 5a^)/24: for neutrino-neutron scattering and 
Cs,p = [4(Cv — 1)^ + 5a2]/24 for neutrino-proton scattering, with Cy = 1/2 + 2sm'^ Ow, 
sin^ 9w ~ 0.23 and a = 1.25 (e.g.; Shapiro & Teukolsky 1983). We calculate the Rosseland mean 
opacity Ks^n, by averaging the neutrino energy over the Fermi-Dirac distribution function with 
zero chemical potential. Using 

^2/ (A;T)2 13.8(A;r)2' ^ ^ 

where k is the gas Boltzmann constant, we find that the total opacity due to neutrino scattering 
is Ks,N = Ks^n + Kg^p = 13.8{Cs,pYp + Cs,nyn)(^op{kT /rrieC^)'^ , where Yn and Yp are the fractions of 
free neutrons n and protons p, respectively. For completely dissociated matter, where this process 
is relevant, the nucleon mass fractions are Yn = 1 — Yf. and Yp = Y^. Here, we take Yg ~ 0.5. This 
leads to a total scattering optical depth given by 



PK'S,N 



H = 2.7 X IQ-^TiVioF, (11) 



'^Note that there is a factor 1/2 difTerence between electron-positron pairs and the neutrino-antineutrino 
contribution, which comes from comparing the electron spin degeneracy (2 spin states) to the neutrino spin degeneracy 
(one helicity state). 
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which has the same temperature and density dependency as equation (8) (when X^^^c = !)• Both 
PWF and NPK found that, for accretion rates above m > 1, where rh = M/(Mqs^^), the disk 
becomes optically thick to its own neutrino emission. This typically occurs within the inner regions 
of the disk. However, neither PWF or NPK studied the effects of including, sclf-consistently in 
the models, the neutrino opacity sources described above. Here, we use a simple prescription to 
include neutrino transport and discuss its effects on models of NDAFs. 



2.3. Basic assumptions and equations 

In this study, (see also NPK), we adopt basic equations based on Newtonian dynamics. As 
in previous work (PWF and NPK), we are interested in the gross properties of the disks and in 
particular on the effects of neutrino transport. For simplicity, we use a steady state disk model. 
Even though accretion in a GRB engine is time dependent, the steady state assumption should be 
a reasonable approximation since the viscous timescales in the inner disk, where all the neutrinos 
processes become important, are much shorter than those in the outer disk, where M is expected 
to vary. 

As in the standard theory of thin accretion disks (e.g., Shakura &; Sunyaev 1973; 
Prank, King Sz Raine 1992), we consider height-averaged quantities and write the 
isothermal sound speed as = P/p and the vertical scale height as H = Cg/^K, where 

= {GM/R'^y/'^ = 2.4 X lO'^mg s~^ is the Keplerian velocity and r is in units of Rs- 

We adopt the standard Shakura-Sunyaev prescription for the kinematic viscosity coefficient, 
v = ac^g/^K: and scale the dimensionless parameter a as a-i = a/O.l (PWF, NPK). The 
continuity equation and an approximate expression for angular momentum balance (e.g., NPK) 
give the following relation for the mass accretion rate: 

M = AirRpHv ^ QirvpH. (12) 

Correspondingly, the radial velocity of the gas is = 2>v j^R. 

In the equation of state we include the contributions from radiation pressure, gas pressure, 
degeneracy pressure and neutrino pressure (see also Popham h Narayan 1995; PWF and NPK). 

\\ . pkT / l + 3Xnuc \ ^ 27r/ic ( 3 ( p\ ^ 

^=12"^ +^l^^j + ^l,w;J UJ +¥' 

where a is the radiation constant and the factor 11/12 includes the contribution of relativistic 
electron positron pairs (as we expect the temperature to be significantly above the pair production 
threshold). In the degeneracy pressure term, /ig is the mass per electron which we take to be 
equal to 2 in agreement with NPK and PWF. Neutrinos also contribute to the equation of state 
and this is taken into account in the forth term, where is the neutrino energy density defined 
as -u^ = l/^aT'^Y^{T^j2 + l/V3)/(r^,/2 + 1/^3 + l/3ra,^,) (Popham & Narayan 1995; see also 
eq.(15)). 
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We write the energy equation as 

, 3GMM _ , , 

where represents the viscous dissipation, q~ = q,J + Q'photo-^ total coohng rate due to 

both neutrino losses {q~) and photodisintegration (gphotoi Eq. 6) and ^adv is the advective cooling 
rate. All the g's correspond to half the disk thickness. For the neutrino flux q~ , we use a simplified 
model for the transport, based on the two-stream approximation and derived by Popham & 
Narayan (1995; see also Hubeny 1990). We write: 

-^y- i7/8)aT^ 

Y(3/4)(r.72 + l/V3 + l/3T.,.J- 

Here t^. = Ta^Ui + Ts^Ui is the sum of absorptive and scattering optical depths calculated for each 
neutrino flavor and Ta^Ui the total absorptive optical depth for each neutrino flavor (e.g. for I'e, the 
absorptive optical depth includes all terms, Ta^^^ = Ta^ueUe + '^a,eN + Ta,brem + Ta,piasmon and the last 
two terms are negligible). Equation (15) gives the correct behavior in the limit of both small and 
large t^^ and Ta,u,^. 

For simplicity, we approximate gadv by, (see e.g., Narayan & Yi 1994; Abramowicz et al. 
1995): 

ds H (n , 3 pk 1+XruA 

where oc (— dlns/dlnr) (see e.g.; Kato, Fukue & Mineshige 1999 for the exact expression for 
^) is taken to be constant and equal to 1. The factor (ll/3)ar^ is the entropy density of the 
radiation. The entropy density of neutrinos, a factor (7/6)ar^ = 4/3 x 7/8aT^, is also added to 
equation (16) when the corresponding pressure term is added to the equation of state. Note that 
degeneracy pressure makes no contribution to the advection term because there is no associated 
entropy. 

We solve numerically equations (13)-(16) to find the temperature and density at a given 
radius and for a given m. We are interested primarily in the properties of the inner regions of 
the accretion flow, where neutrino processes are important. We thus concentrate on the range of 
radii from rmin = 3 to rmax = 200. (From the work of PWF and NPK we know that, for r ^ 100, 
the flows are fully advection dominated since neutrino cooling is not important and photons are 
completely trapped). 



^Equation (15) was derived in the context of radiative transport and assumes that opacities and emissivities are 
independent of the z coordinate. Although we do not solve for the vertical structure of the flow, this assumption 
is less accurate for the neutrino transport where all cross sections are a function of temperature and hence of the 
vertical disk structure 
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3. Numerical results 

3.1. Gas profiles, neutrino optical depths and advection 

Temperature and density profiles calculated from our model are shown in Figure 1 and the 
corresponding pressure components in Figure 2. We show our solutions for three values of the 
accretion rate, m = 0.1, 1, 10 (long dashed, solid and short dashed line, respectively). The thinner 
lines show the corresponding PWF solutions (see also their Figure 3). For easier comparison 
between the two sets of solutions, we extend our solution range out to r = 1000. PWF and NPK 
have shown that, at large radii, densities and temperatures are too small for neutrino cooling to 
be significant, while optical depths are too large for photon cooling so that the flows are simply 
advection dominated. Therefore the solutions should not be sensitive to any assumption made 
about neutrino cooling and should compare well in this region. Figure 1 shows that, at large 
radii, there is indeed good agreement, both in the slope and in the normalization, between our 
solutions and those of PWF. In this region, radiation pressure dominates for rh < 10, whereas 
the degeneracy and gas pressure components become increasingly more dominant for the higher 
values of m (Figure 2). In the region from r ~ 100 to ~ 300, the temperature profiles flatten. 
Here, photodisintegration and, more importantly, neutrino emission starts cooling the gas. The 
treatment of neutrino transfer inside this region becomes important and the two solutions start 
differing. Note that the onset of photodisintegration produces a sharp feature in the PWF solution 
in this range of radii. Such a feature is also reproduced in our calculation although it is somewhat 
weaker, particularly as seen in the density profiles. 

Inside radii of r ~ 30 — 100, and for m > 1, our temperature solution steepens again whereas 
the PWF profile flattens. A similar trend is observed in the density profiles. In this region gas 
pressure becomes the dominant term in the equation of state (Figure 2). 

Figure 3 shows contour plots for three quantities: the optical depths Ta^i, = J2'''a,i/i 
(dotted lines) and Tg^u = J2 '^s,vi (dashed Unc) in the top panel, the advection 
parameter, / = qadv/o^ = 1 ^ <? Z*?^ in the middle panel and the viscous timescale, 
tacc = /i^ = (H / R)^'^Q]^^ in tlic bottom panel. These quantities have been calculated for 
m ranging from 0.1 to 10 in a region of the disk from r = 3 to 200 (where neutrino cooling is 
important) . The top panel of Figure 3 shows that the optical depth is dominated by the absorptive 
opacity in the inner regions of the flow (and in particular by the inverse of proton and neutron 
capture, given by Eq. 4). The optical depth exceeds 1 for r in the range ^5 — 30 and m = 0.1 — 10. 
Note that the cooling function in equation (15) reduces to the optically thin expression for small 
optical depths but it differs significantly from the latter at optical depths ~ 1 and larger. The 
treatment of neutrino transport via equation (15), in the optically thick regions, is what gives rise 
to the major differences between our work and that of PWF. 

The most significant consequence of the gas becoming opaque is that neutrinos are trapped in 
the inner regions of the flow and the neutrino emission is partially suppressed. This is illustrated 
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by the middle panel in Figure 3. Energy advection becomes important (i.e. / > 0.5: advection 
dominates over cooling) in the inner region of the flows and in particular in the region where the 
optical depths become significantly larger than 1. By comparing the contours of the optical depths 
and / one finds that the region where most of the neutrino cooling occurs (where / decreases) 
lies roughly along (or just outside) the r = 1 contours. With increasing rh the neutrino cooling 
becomes increasingly less efficient as the r = 1 surface, moves further out in the flow. Consistently 
with previous work, we also find that / (i.e. advection) also dominates at large radii. In particular, 
/ ~ 1 for r > 100; this is equivalent to the statement that the cooling timescale is much longer 
than the accretion timescales, so the energy is advected inward before it can be radiated away. 

In our solution, which is mostly advection-dominated, we have H ^ R. Therefore, tacc is 
almost independent of m and approaches the simple dependence: tacc ~ l/a^K- Our values of 
tacc are typically a factor ~ 10 smaller than in Figure 1 of NPK, where the advection term and the 
appropriate opacities were not taken into account self-consistently within the NDAF region. 

In the optically thin approximation used by PWF, the neutrino emission increases with 
increasing m and energy advection is not important in the inner regions of the flow. Although 
PWF also estimated that the flows become optically-thick to their neutrino emission above m ~ 1 
(consistent with our results), they did not account for neutrino opacities in their models. Our 
results, although based on a simple model, emphasize the importance of accounting for neutrino 
transfer^. We find that, for most of the range of accretion rates of interest for GRB models, a 
neutrino-dominated accretion flow becomes optically thick to neutrinos, and advection of energy 
provides the most significant cooling term in its central regions. In §5, we will discuss the 
implications of these results for the accretion flow neutrino luminosities in GRBs. 



NPK discussed the stability properties of their NDAF solution. Since our model differs 
considerably from theirs, we repeat the analysis here. 

The general condition for thermal stability is: 



where Q = q + Qadv The top panel of Figure 4 shows the two functions and Q (solid 
and dashed line respectively) as a function of the temperature of the gas. The radius is fixed at 



^We note that PWF considered a fuUy-relativistic Kerr geometry. Our solution does not take into account the 
effects of black hole rotation. However the differences between our solution and that of PWF are always much larger 
than those introduced by allowing for black hole rotation (see e.g.; Figure 6 in PWF). Although the black hole is most 
likely rotating, we do not consider our main conclusions to be significantly limited by adopting the simple Newtonian 
treatment. This is because, as shown in Figure 3, most of the neutrino energy comes from radii of several Rs- 



4. Stability 




(17) 
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r = 100 (for the lower set of curves) and r = 10 (for the upper set of curves which have also been 
renormalized by a factor 10^), while the surface density is taken to be E = lO^^gcm"^. The 
equilibrium value of T is given by the intersection of the two curves. We see that the stability 
condition above is satisfied at this point, so that the solution is thermally stable. 

The condition for viscous stability is given by 

i > <^«) 

The bottom panel of Figure 4 shows m versus S at r = 100 (solid line) and at r = 10 (dotted 
line). This shows that the disk is viscously stable. We have checked both the thermal and viscous 
stability conditions for a variety of r and m and we find that the fiows are stable under all 
conditions. This is a well known property of flows in which advection is important. 

Finally, we also check for the gravitational stability condition, namely that the Toomre 
parameter Qt should be larger than 1. For a Keplerian disk, Qt is given by: 

where n is the epicyclic frequency. Figure 5 shows the variation of Qt with radius for 
m = 0.1, 1, 3, 10. We see that Qt decreases with increasing r, so that the flows are most unstable 
on the outside (as found also by NPK). However, only for the largest m ~ 10 and for r ^ 50 does 
Qt go below unity, signifying gravitational instability. Thus the flows are gravitationally stable 
under almost all conditions of interest. 



5. Discussion 

The relevance of accretion processes in the central engine of GRBs was highlighted by 
Narayan et al. (1992). PWF carried out the first detailed study of accretion around black holes 
with ultra-high m (of order a solar mass per second), and identified a new mode of accretion which 
they named neutrino-dominated accretion flow (NDAF). In a subsequent study, NPK noted that 
accretion flows with such high m are highly advection dominated at large radii (r ^ 100). They 
suggested that the flow at these outer radii may take the form of a convection-dominated accretion 
flow (CDAF; Narayan, Igumenshchev & Abramowicz 2000; Quataert & Gruzinov 2000) or a 
related kind of flow (e.g., Blandford &: Begelnian 1999). As a result, only a small fraction of the 
available mass accretes on the black hole, the bulk of the gas being ejected from the system. Close 
to the black hole (r ^ 100), however, NPK confirmed PWF's result, that a cooling-dominated 
NDAF should be present. 

PWF and NPK made the simple assumption that the accreting gas is optically thin to its 
own neutrino emission. We have improved on this by including the effects of neutrino transfer 
(via a simple prescription) and the neutrino opacities self-consistently in the model. With these 



- 11 - 



improvements, we find that the central regions of NDAFs are typically opaque, so that the 

neutrinos emitted by the accreting gas are largely trapped. The optically-thick region extends out 
tor~4 — 5atm~ 0.1 to r 30 — 40 at m 10 (Figure 3). We have shown that, above rh > 1, 
the majority of the energy liberated by viscous dissipation is advected with the flow instead of 
being emitted in the form of neutrinos, and energy advection is much more important in the inner 
regions of the accretion flows than previously realized (Figure 3) . 

In order to further explore the implications of these results we now calculate the neutrino 
luminosity from our model and the fraction of the binding energy carried out to infinity by 
neutrinos. We also derive the neutrino Eddington luminosity for such flows. 

Figure 6 shows the neutrino luminosity from the accretion flow, = f^^.^ 2'Kq~rdr, (in units 
of 10^^ ergs""^ ; solid line; upper panel). Also, using the opacities discussed in §2.2 we derive the 
neutrino Eddington luminosity of the flow: 

Lmd,u = ^^^^^ ~ 9 X lO^^^msTfi^ erg s"! , (20) 

where, for illustration the numerical value given on the right uses only the two dominant 
components of the opacity: neutron and proton absorption and scattering (for electron neutrinos 
this gives = Kue,a + K'Ue,s = 5.5 x 10~^^rfi cm^ g^^; Equations 8 and 11; note that for the 
calculation shown in Figure 6 all opacity terms are used). Because the neutrino cross sections are 
a function of neutrino energies, the neutrino Eddington luminosity is a function of both black hole 
mass and gas temperature. In Figure 6, we show the neutrino Eddington luminosity calculated for 
Til = Tii{r = 3), the temperature at the inner radius of the flow; at larger radii the Eddington 
luminosity increases. In the middle panel of Figure 6 we also show the neutrino radiative efficiency, 
defined, as in standard accretion theory, by rj^, = L^/Mc?. 

The neutrino luminosity from the flow is ~ 10^^ ergs~^ at m ~ 0.01 and increases almost 
linearly with m up to m ~ 0.1 (Figure 6). Between 0.1 < m < 1, the luminosity flattens off 
significantly, ranging within 6 - 8 x 10^^ erg . Above m ~ 1 the neutrino luminosity stays 
virtually constant at ~ 8 x lO^^ergs"^ . In contrast, in the PWF solution, which assumes 
optically thin neutrino emission, the neutrino luminosity ranges from ~ lO^'^"^^ ergs~^ for 
m = 0.1 to up to Ly 10"'''* crgs~^ for m = 10 (see their Table 3). Thus, one important 
consequence of the addition of neutrino opacities is that the flow luminosity varies only by a factor 
of a few over the range of accretion rates appropriate for popular models of GRB progenitors. The 
luminosity of the flow is almost constant because the efficiency with which energy is transported 
out of the flow by neutrinos, r/j, (bottom panel. Figure 6) decreases with increasing m (or 
equivalently, as shown in Figure 5, decreases as neutrinos become trapped and energy advection 
becomes more dominant at larger m). We show that t] 0.1 (consistent with the efficiency 
expected from a geometrically thin, cooling dominated Newtonian accretion disk) up to m ~ 0.1 
and decreases (almost linearly) with increasing m (reaching r] ~a few 10-3 at m = 10). The 
fraction of energy transported away by neutrinos decreases as the inner regions of the flow become 
neutrino opaque (see Figure 3). 



- 12 - 



Figure 6 also shows the neutrino Eddington luminosity evaluated at the inner edge of the disk 
(dashed line). In the inner regions of the flow, the temperature is the highest, and consequently the 
Eddington luminosity the lowest (Eq. 20), so that ~ L-^Ad^v for m > 1. Momentum deposition 
in the inner regions of the flow by neutrinos may therefore be quite effective as a mechanism 
responsible for ejection. 

We now examine the implications of our results for GRBs. A number of authors, including 
Eichler et al. (1989; see also Narayan et al. 1992) suggested that vu annihilation around merging 
compact binaries might produce a relativistic e^e~ fireball with sufficient energy to power a GRB. 
Detailed numerical simulations by Janka et al. (1999) show that this process might conceivably 
power short GRBs (those with durations under a couple of seconds or so), provided that there is a 
modest level of beaming. However, we have shown that neutrino advection is important in most 
of the parameter space of these models, which may have a significant impact on the efficiency of 
this process. We estimate the luminosity due to vi> annihilation along z-axis above the disk (using 
an improved version of Eq. 10 in Ruffert et al. (1997) provided by Thomas H. Janka, private 
communication) to be: 

L,p ~ 6 X 10-=^^—- -2 -— / de ^ (cos - cos e,^,if erg s'^ , (21) 

where the constant in front takes into account the neutrino cross sections. Li, is given above 
and shown in Figure 6 and i?surf and (Ei^) are evaluated at the neutrinosphere, defined as the 
surface at which the emergent neutrinos originate. Hence, {E^) is the average energy of the 
escaping neutrinos given by (E^) = 3.7 kTu ^surf (Eq. 10) with Tn^surf = 7'ii(?'surf = r{Tu = 2/3)) 
and Tu = Ta,u + Ts,u- If the spectrum is a blackbody, the radius of this surface is located near 
and above the layer of optical depth = 2/3 (this is because the opacity is mostly dominated 
by absorption processes). The contour of t^, = 2/3 and the radius of the neutrinosphere 
as a function of m is shown in Figure 3. We find that the radius of the neutrinosphere 
to be given by, rgm-f = r^=2/3 ~ 17m^/^ in the range 0.1 < m < 10. At this radius the 
temperature of the neutrinosphere is roughly Tn^surf ~ 0.2m~'^/'^^ which is fairly insensitive to 
rh in the range from 0.1 ^ m < 10). The integral part in equation (21) is a geometrical factor, 
where $ = 3/4[l - 2{fi)^ + (^2)2 ^ i^2{l - (/i^))^], where (^") = i^^", d^l^l'^ / ^^^^^^ d^l with 

^surf = cos6lsurf = l/\/l+7^irf/^> /^min = cosO^^^ = l/-y/l + {rl^^J z'^) and e = z/rsurf. 

The bottom panel of Figure 6 shows the neutrino annihilation luminosity, Ly^ in units of 
10^^ erg s ^ . Li/jy increases up to its maximum value of ^ 10^^ erg s ^ at tti ^ 1 and slightly 
decreases for larger rh. Our estimate of Ly^p at m = 1 agrees well with values estimated in the 
simulations of Ruffert & Janka (1999) and PWF (although our results are less accurate than those 
of PWF for m < 0.1). In accordance with PWF we find that energetic events can only be achieved 



^"To derive the analytical scalings we solve Eqs. (13) and (14) assuming P — Pgas and g"*" = q^^^- We have shown 
that is a good approximation within the optically thick region. 
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for rh > 0.1 (below which the neutrino annihilation efficiency decreases very sharply; see also 
their Table 3); but in contrast with their results we do not find that increasingly more energetic 
events can be achieved for larger accretion rates. Our calculations imply that the efficiency of iyi> 
annihilation remains constant (or slightly decreases) for m > 1. This is expected, as we have found 
that neutrinos are increasingly more trapped in the disk. (Our estimate may be uncertain by a 
factor ~ 2; this is mostly due to the fact that our solution gives only height-averaged quantities 
whereas the vertical disk stratification may be important in this calculation). Note also that recent 
results from hydro dynamical calculations carried out by Lcc &; Ramirez-Ruiz (2002) show that 
annihilation can only produce bursts from impulsive energy inputs, as the annihilation luminosity 
scales as This further restricts the importance of this process to a small fraction of bursts. 

Energy extraction from the disk is also possible by MHD processes. Such processes are 
broadly based on the expectation that the differential rotation of the disk will amplify pre-existing 
magnetic fields, until they approach equipartition with the gas kinetic energy. Proposed 
mechanisms include Parker instabilities in the disk leading to reconnection, relativistic flares 
and winds (Narayan et al. 1992; Meszaros & Rees 1997) or the Blandford-Znajek mechanism 
(Blandford & Znajek 1977; hereafter BZ). Of these, perhaps, the Blandford-Znajek efficiency is 
the easiest to estimate (at least roughly). Wc follow the common assumption that the magnetic 
field in the disk will rise to some fraction of its equipartition value B'^/Stt ~ pc^. Typical values of 
pel are 10^°~^^ ergcm"^, for 0.1 < m < 10 (see Figure 2), implying a field strength of ~ 10^^~^^ G 
if we make the conservative assumption that B is only 10% of its equipartition value. Consistent 
with earlier work (PWF; RufFert &; Janka 1999), the BZ jet luminosity at m = 10 is then 

L..^(^!^y^iA>^,^V(^)\^f e.,.-, (22) 

where Rh = 2GM/(? is the black hole radius and a = Rh^h/c is the dimensionless black hole spin 
parameter (0 < o < 1). Lbz about 2 orders of magnitude larger than the neutrino annihilation 
m ~ 10. For larger equipartition fractions, the ratio Lbz/Lup becomes large even at 
m ~ 0.1 or 1. It is obvious that the energy liberated by this mechanism is simply proportional 
to the internal energy density of the disk and hence is not affected by advection. Indeed, Livio, 
Ogilvie & Pringle (1999) have argued that the BZ mechanism might be most relevant for advection 
dominated flows (see also Meier 2001). Numerical simulations such as those by MacFadyen & 
Woosley (1999) and MacFadyen, Woosley k, Heger (2001) have shown that the BZ mechanism 
would be an efficient mechanism capable of liberating a large fraction of the black hole spin energy. 

By similar electromagnetic considerations as those used to derive the BZ luminosity, Livio et 
al.(1999) estimate that the electromagnetic/ wind power output from a disk is given by: 

where, following Livio et al., we take R^^d/c ~ (Rh/Rd)^^^- -Bpd is the poloidal field in the disk 
given, approximately, by Bpd ~ {H/R)B (see also Merloni & Fabian 2002). Because H/R ~ 1 in 
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the inner region of the disk we hayeBpa B. Rd is a factor of a few times Rh (above we take 
R(i as the radius of the innermost stable orbit, at which the disk energy density is the highest). 
Since a? < 1 and, Rd is at least a few times Rh, it is evident that in a geometrically thick disk, 
the electromagnetic/wind luminosity from the disk can easily be larger than the BZ luminosity, 
hence Ld » Lyp- Our calculations therefore indicate that, with increasing accretion rates, 
MHD processes become significantly more efficient at releasing energy than neutrino annihilation 
processes, hence they are probably the most viable mechanisms for energy extraction in these 
systems. 

Although the stability analysis presented in §4 shows that the accretion flows we have studied 
are intrinsically stable, we note that the viscous time scale (Fig. 3) can be as short as a small 
fraction of a second. Therefore, if the accretion disk is fed in a variable manner, e.g., via fallback 
material after a supernova explosion as in the collapsar model, we may expect variations in the 
accretion rate. This might explain the complex light curves of GRB in at least some of the 
scenarios discussed above (e.g., MacFadyen & Woosley, 1999). 
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Fig.l — Temperature and density profiles in the top and bottom panel respectively. The 
profiles are shown for three values of the accretion rate: m = 10 - short dashed lines; m = I - solid 
lines; m = 0.1 - long dashed lines. The PWF numerical solutions (thinner lines - see also their Fig. 
3) are also plotted for direct comparison. 

Fig. 2 — Pressure components for three values of rh. The gas pressure is shown by the solid 
line, degeneracy pressure by the dotted line, radiation pressure by the dashed line, neutrino 
pressure by the long dashed line 

Fig. 3 — Top panel: Contours showing the values of scattering optical depth (r^ - dashed 
line) and absorptive optical depth (ra - dotted line) for a range of m = 1 — 10 and r = 3 — 200. 
The T = 2/3 contour is shown with the solid line. Middle panel: Contours showing the advection 
parameter / = Qadv/Q^ for the same range of m and as a function of r. For most regions of 
parameter space / > 0.5 and advection dominates. Bottom panel: contours of the accretion time, 

tacc- 

Fig.4 — Stability analysis. The top panel shows (dashed line) and Q~ = q~ + Qadv (solid 
line) versus T for r = 100 (bottom curves) and r = 10 (top curves - renormalized by a factor of 
10^ for better clarity) and for S = 10^^ g cm~^ . According to the criterion given in Equation 10, 
the flow is thermally stable. The lower panel shows rh versus S for r = 100 (solid line) and r = 10 
(dotted line). The viscous stability condition (Eq. 11) is always satisfied. 

Fig. 5 — The Toomre parameter Qt (Eq. 12) as a function of radius for the four values of the 
accretion rate as in Fig. 2. Qt is always > 1 (for m > 3) implying that the flow is gravitationally 
stable throughout and unstable in the outer regions for m > 3. 

Fig.6 — Upper panel: the solid line represents the neutrino luminosity and the dashed lines 
the neutrino Eddington luminosity as a function of the accretion rate and for two values of the 
effective temperature. Middle panel: the solid line shows the radiative efficiency of the accretion 
flows as a function of m. Bottom panel: Li,p as a function of rh. The dotted line below rh = 0.1 
represents the most uncertain (upper limit) estimates. Better values of Li,p for m < 0.1 are given 
in Table 3 of PWF. 
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